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ABSTRACT 

Planet formation is profoundly impacted by the properties of protoplanetary 
disks and their central star. However, how disk properties vary with stellar 
parameters remain poorly known. Here we present the first comprehensive, com- 
parative Spitzer/IRS study of the dust and gas properties of disks around young 
sun-like stars (K1-M5) and cool stars/brown dwarfs (M5-M9). The compari- 
son of these two large samples of over 60 sources reveal major differences in the 
evolution of both the dust and gas components. 

We report the first detection of organic molecules in disks around brown 
dwarfs. The detection rate statistics and the line flux ratios of HCN and C 2 H 2 
show a striking difference between the two samples, demonstrating a significant 
under-abundance of HCN relative to C 2 H 2 in the disk surface of cool stars. We 
propose this to originate from the large difference in the UV-irradiation around 
the two types of sources. The statistical comparison of the 10 ^m silicate emission 
features also reveals a difference between the two samples. Cool stars and brown 
dwarfs show weaker features arising from more processed silicate grains in the 
disk atmosphere. 

These findings complement previous indications of flatter disk structures and 
longer disk lifetimes around cool stars. Our results highlight important differ- 
ences in the chemical and physical evolution of protoplanetary disks as function of 
stellar mass, temperature, and radiation field which should be taken into account 
in planet formation models. We note that the different chemistry of pre-planetary 
materials in the disk may also influence the bulk composition and volatile con- 
tent of the forming planets. In particular, if exogenous HCN has played a key 
role in the synthesis of prebiotic molecules on Earth as proposed, then prebiotic 
chemistry may unfold differently on planets around cool stars. 
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Subject headings: 
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Introduction 



Circumstellar disks are a natural outcome of the star formation process and the 
birthplace of planetary systems. Their evolution and their lifetime determine what type 
of planets can form. Most studies so fa r have focused on characterizing protoplanetary 
disks around young sun-like stars (e.g. iMeyer et al. 



2007 for a review). It is now well 



established that most sun-like stars clear out their primordial dust disk in 5Myr and only 



Hernandez et al. 



20071 ). The 



a few percent of them still retain it by an age of 10 Myr (e.g. 
dispe rsal of the gas component is less well character ized but appears to proceed equally 



fast (IPascucci et al. 



2006 



Sicilia-Aguilar et al 



20061 ). This rapid clearing of dust and gas 



around sun-like stars matches the fast growth of planetesimals and gas dispersal in the 



protosolar nebula as inferred from meteorites, asteroids, and planets (IPascucci fc Tachibana 



20091 ). But do disks around stars of different masses evolve similarly? Recent studies hint 



that disks arou nd cool/late-type M stars, the most typical protoplanetary disks, evolve 



differently (e.g. 



Apai et al. 



20081 ). There are at least three major differences. 



F irst, the dispersa 



stars. ICarpenter et al. 



of p rimordial dust disks seems to operate less efficiently for cool 
(120061 ) find that the disk frequency of low-mass stars (~ 0.1 — 1.2 M ) 



is statistically higher than that of high er-mass stars i n the 5 Myr-old Upper Sc o OB 



association. I nfrare d excess studies by 



Sterzik et al. 



(120041 ) 



Scholz et al. 



(120071 ). and 



Riaz &: Gizisl (120081 ) all suggest that this trend extends to the brown dwarf regime. Second, 



there is evidence that disks around cool stars ha ye, on average, flatter disk structures 



than the more fla red disks around sun-like stars (IPascucci et al 



Allers et al. 



2003 



Mohanty et al. 



2004 



Apai et al. 



2004 



Bouy et al 



20081 ). 



20061 but see also the flared disks from 
Third, the dust probed through the 10 fim silicate emission feature appears to be more 
processed around cool stars tha n around sun-like stars of similar age, sugge sting more rapid 



grain growth in cool star disks (lApai et al. 



2005 



Kessler-Silacci et al. 



20061 ). 
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However, there has not been a detailed comparison of the dust properties of disks 
around sun-like and cool stars and no investigation of their gas content. Here we present 
the first such comparative study. We will show that there are statistically significant 
differences between cool and sun-like stars in their gas and dust disk properties. First, we 
will demonstrate that cool stars/brown dwarfs have much weaker continuum-subtracted 
10 /jiii features than higher-mass stars (Sect. 13. ip implying different dust populations at the 
observed disk radii (Sect. 14. ip . Second, we will show that HCN emission is often present in 
sun-like star spectra, but absent from the cool star spectra. Similarly, the median sun-like 
star spectrum has an HCN/C2H2 flux ratio that is an order of magnitude higher than the 
median cool star spectrum (Sect. BHD - We will see that this arises from an under-abundance 
of HCN in the disk atmosphere of cool stars (Sect. H~2l) . These findings highlight how stellar 
mass, and radiation field affect the physical and chemical evolution of protoplanetary disks 
including the organic compounds available during planet formation. 



2. Targets, Observations, and Data Reduction 



2.1. Target Selection 



We compare two large samples of disks around sun-like and cool stars (including 
brown dwarfs) that have Spitzer spectra obtained with the same resolution and very similar 
signal-to-noise. The Taurus and the Cha I star-forming regions provide the best samples 
of protoplanetary disks to carry out our comparative study. These low-mass star forming 
regions are at about the same distance and do not contain early B or O type stars. In 



add ition, their median age (~1 and 2 Myr' 



e.g. 



Kenyon et al. 



2008 



and 



Luhman et al. 



is the same within our ability to discern (see, 



2008bl for a review on Taurus and Cha I) . 



The sun-like star sample is the same as in 



Pascucci et al. 



(120081 ): it contains 23 single 
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and 21 binary stars from the Taurus- Auriga star-forming region with spectral types between 
K1-M5, temperatures betw een 3,000-5,000 K, and stellar masses between 0.4-2 M (see also 



Kenyon fc Hartmann 



19951 for other stellar properties). Two-thirds of the binary systems 



have stellar companions between 0.1"-1", with a mean projected separation of 0.4", or 56 AU 



at the distance of the Taurus star-forming region. The low-resolution IRS spectra 
sun-like star sample were acquired as part of the IRS/ GTO prog ram (IFurlan et al. 



and have been reduced as described in 



Pascucci et al. 



Pascucci et al 



(120081 ) and 



Bouwman et al. 



or the 



2006) 



20081 ) 



(120081 ) showed that there is no statistically significant difference between 
the strength of the 10 /mi emission features from the single and binary systems in these 
two samples. Similarly, they showed that the ensemble of spectral energy distributions is 
indistinguishable between the single and binary system samples. 

Our cool star sample includes objects with spectral types between M5-M9, stellar 



temperat ures between ~2 , 500-3, 00C 



and also 



Apai et al 



2005 



Luhman 



< an d masses between ~0. 04-0.1 5 M (see Tabled] 



20071 ). These objects were selected from the Cha I 



star-forming region following two main criteria: i) reliable spectral clas sification, and ii ) 



mid-i n frared excess emission. W e used the optic al/infr ared surveys b y 



(I2000f ) 



Lopez-Marti et al. 



(J200J), iLuhmanl (120041 ) . and 



Luhman et al. 



Comeron et al. 



((2005J) to identify 



sources close to or below the substellar boundary. We then correlate d these sources with the 



ISO/ISOCAM de tections at 6.7 and 14 3 /m from 



8 /mi fluxes from 



Luhman et al. 



Per si et al. 



(I2000f ) and the Spitzer/IRAC 



(l2008al ). The ISOCAM-selected sample consists of eight 
sources: Hn 2, CHXR 15, Cha Ha 1, Cha Ha 9, Cha Ha 2, ISO 138, Cha Ha 6, and 
ISO 217. These sources were observed in March 2005 as part of our Cycle- 1 Spitzer proposal 
and were selected to have fluxes larger than 10 mJy in at least one of the two wavelengths 
to ensure a reliable detection of excess emission. Another eleven sources were selected later 



based on excess emission in the IRAC 5.8 and/or bandpasses (iLuhman et al 



2008a 



and were observed as part of our Cycle-3 Spitzer proposal in August 2006 and March 2007. 
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The IRAC-selected sources are: ESO Ha 559, ISO 79, CHSM 9484, Cha 23943, ISO 147, 
Cha 9086, T37, ISO 165, ISO 252, Hn 13, and Cha 21330. 



2.2. Cool Star Sample: Observations 



We used the infrared spectrograph (IRS, iHouck et al. 



20041 ) on-board the Spitzer 



Space Telescope to perform low-resolution (R ~64-128) spectroscopy of the 10 jttm silicate 
emission feature of the selected cool stars. We employed for all sources the SL1 module 
covering the wavelength range between 7.4-14.5 /im. For several sources of the Cycle-1 
proposal we also adopted the LL1 module (19.5 - 38.0 /xm) aiming to probe solid-state 
features at long wavelengths. Ho wever, because m ost of the selected brown dwarfs turned 



out to have flat disk structures (lApai et al. 



20051 ). the signal-to-noise ratio of the LL1 



observations is only modest. This means that we can use the LL1 data only to trace 
the infrared continuum. In view of this finding, we preferred to drop the LL1 module 
for the Cycle-3 observations. We chose instead the SL2 module (5.2-7.7 /im) to better 
establish the continuum slope around the 10 /im silicate emission feature. The parameters 
of the observations are summarized in Table [3 Most of the sources were centered in the 
spectrograph slit using the high-accuracy PCRS mode (1 a positional uncertainty of 0.4"). 
The exceptions are Cha 9086, ISO 217, and Hn 13 where no suitable optical reference star 
was available. For these sources, we used the medium accuracy infrared peak-up array 
integrated onto the IRS focal plate providing a 1 a positional uncertainty of 1". All peak-up 
operations were successful and no significant flux loss is expected from any of our sources. 
Each source was checked individually to confirm that only one point source was covered 
by the slit. We used fixed cluster offset-mode observations and grouped nearby targets 
to reduce the telescope overhead. By using the staring mode observation template, we 
acquired each target in two nod positions along the spatial direction of the slit, at 1/3 and 
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2/3 of the slit length. 



2.3. Data Reduction 



The sun-like and cool star samples have been reduced in an identical way, with only 
ve ry minor differences. We have already published the spectra of the sun-like star sample 



in 



Pascucci et al. 



(120081 ). Here we adopted the already reduced spectra from our previous 



paper and will now only discuss the reduction of the cool star sample. 

The cool star data have been processed through the SSC pipeline version S15.3.0 and 
S16.1.0. Our goals required a data quality higher than that produced by the SSC pipeline. 
We start our data reduction from the so-called droopres products and use the SMART 



reduction package (IHigdon et al. 



2004 ). in combination with IDL routines developed for 



the FEPS Spitzer Scienc e Legacy program 



are outlined in detail in 



Bouwman et al 



( Meyer et al. 



series of publication s (see, e.g. iMeyer et al 



20061 ). The data reduction steps 



(2008) and have been successfully used in a 



2004 



Pascucci et al. 



Apai et al. 



2005 



Bouwman et al. 



2006 



20081 ). In brief, we first subtracted the pairs of imaged spectra acquired 



along the spatial direction of the slit to correct for background emission, stray light, and 
pixels with anomalous dark current. We replaced pixels flagged as bad by the SSC pipeline 
(and additional 3-4 pixels per frame identified by eye) by interpolating over neighboring 
good pixels. Spectra were extracted from the background-subtracted pixel-corrected images 
using a 6-pixel fixed-width aperture in the spatial direction. 

Because all targets have been acquired with a high or moderate peak-up option 
(pointing accuracy within 0.4-1") and the majority are too faint to perform source profile 
fitting, we opted to fix the position of the aperture for each spectral order to that derived by 
the FEPS legacy program. After extracting the spectra for each order, nod, and cycle, we 
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computed a mean spectrum for each order and as uncertainty we quoted the la standard 
deviation of the distribution of the flux densities measured at the given wavelength. We 
converted our output sp ectra to flux by applyi ng; the spectral response function derived by 



the FEPS Legacy team (jBouwman et al. 
the quoted uncertainties. 



20081 ). and propagate the calibration error into 



3. Results 



The SL1 observations detected all cool stars and achieved a sub-mJy accuracy. 
We confirm the ISOCAM excess emission detections in all but two sources: Our 
Spitzer/IRS observations of Hn 2 and CHXR 15 show flux densities fully consistent with 
photospheric emiss i on an d we identify no hints of the ISOCAM excess emission reported by 



Lopez-Martf et al. 



(120041 ) . It seems likely that these two sources have been missclassified: 



photospheric emission from sources earlier than our spec tral type limit of M6 can fully 



explain the ISOCAM observations (jLuhman et al. 



2008al . see also Table [T]) . We excluded 



these two sources from the following analysis and discussion. In addition, we find that 



ISO 1 38 has IRS excess emiss ion inconsistent with the ISOCAM photometry ( 



2000 



Lopez-Martf et al. 



20041 ) , but fully consistent with the IRAC 8 /xm by 



Persi et al 



Luhman et al. 



(j2008al ). see also Fig. [2j The IRS and the IRAC fluxes at 8 fim are about eight times lower 
than the ISOCAM flux at 6.7 ^an. We suspect that the source 2 MASS 11082238-7730277 a t 
~18" from ISO 138 has contaminated the ISOCAM estimates by 



Lopez-Martf et al. 



(120041 ) . 



The LL1 observations detected all disks in our sample. The low LL1 fluxes confirm the 



previ ous suspicion that most very low mass objects have flatter disk structures (lApai et al. 
20051 ). However, given the low S/N these observations do not further our understanding of 
the dust composition and will not be discussed in the rest of the paper. 
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We find that 14 out of 17 cool star disks have 10 /im silicate emission features. These, 
combined with the sun-like stars and the Herbig Ae/Be stars as described in Sect. 13.11 
make a sample of 72 disks where the dust properties can be studied as a function of stellar 
mass/lumin osity. To compute the strength of the 10/im features we followed the procedure 



described in 



Pascuccietal 



( 120081 ). The continuum below the feature is estimated as follows. 
For the sources that have IRS spectra covering the wavelength region between 5.3-14.2 /im 
we fit a third-order polynomial between 6 and 8 /im and between 12 and 14 /mi. These 
sources comprise all the sun-like stars and the following cool stars: ESO-Ha559, ISO 79, 
CHSM 9484, ISO 147, T37, ISO 165, ISO 252, and Hn 13. For sources that have a smaller 
wavelength coverage we fit a first order polynomial between the minimum wavelength and 
8.2 //m and between 12 /im and the maximum wavelength. The sources belonging to the last 
sample are all the Herbig Ae/Be disks and the following cool stars: Cha Hal, Cha Ha9, 
Cha Ha2, ISO 138, Cha Ha6, and ISO 217. 

We report the detection of the v 5 vibration-rotation band of acetylene (C2H 2 ) around 
13.7 /im and of the v-i band of hydrogen cyanide (HCN) around 14.0 /im toward at least 20 
disks in the combined sample of sun-like and cool stars (Sect. 13.21) . We identify the v 5 and 
v-i bands of C2H2 and HCN using the following procedure. First, we fit a linear continuum 
in a wavelength region outside the expected emission bands (between 13.4-13.60/im and 
14. 15-14. 3/im). Then, we integrate the flux of the continuum-subtracted spectra within the 
wavelength region of the low-resolution spectrum where the emission bands are expected to 
be (13.6-13.81/mi for C2H2 and 13.81-14. 15/im for HCN). The uncertainty on the measured 
band flux is dominated by two errors: a) the error on the continuum subtraction, which 
we compute as the standard deviation of the pixels outside the emission feature in the 
continuum-subtracted spectrum; and, b) the propagation of the measurement errors at each 
pixel. We classify emission bands as "firm" detections or "possible" detections, based on 
the S/N values derived above. We successfully confirmed this classification by individually 
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inspecting every spectra by eye. In the case of the cool star sample possible detections have 
S/N between 2 and 5 while all firm detections have a S/N > 7. In the case of the sun-like 
sample the S/N estimate is more uncertain for two reasons: i) the fit of the continuum is 
less certain because the continuum is less flat than that from cool star spectra; ii) the errors 
at each wavelength are more uncertain because the sun-like stars have only 4 exposures per 
nod position (for comparison cool stars have more than 10 exposures per nod). Most firm 
identifications have S/N greater than 3 in the sun-like star sample (possible identifications 
have S/N slightly lower than 3). Five out of 17 cool star disks have "firm" detections 
of C2H2, other 4 sources have possible detections. Only 1 cool star shows a "possible" 
detection of the HCN rovibrational emission band. Among the 44 sun-like stars, there are 4 
"firm" and 3 more "possible" detections of C2H2, while there are 13 "firm" detections and 
7 more "possible" detections of HCN (see also Table [3]). 



3.1. Weak Silicate Emission Features in Cool Star Disks 

Grains from the interstellar medium (ISM) provide the initial solids that comprise 
forming circumstellar disks. These grains are likely a mix ture of silica t es and carbon with a 
range in size from a few tens of A up to ~0.2 /an (see, e.g. 



Draine 



2003 



Gail fc Hoppe 



20091 ). 



The 10 /im absorption features toward several lines of sight show that mos t of the silcate s 



are a morphous and crystalline material constitutes less than 2% by mass (IKemper et al. 



20051 ). These primordial grains undergo major chemical and physical changes during and 
after the formation of circumstellar disks, some of which are reflected in the 10 /iin silicate 
emission feature. Silicate emission features from our targets show a considerable diversity 
ranging from the strong, sharp feature of Cha Ha 1 to the almost featureless spectra of 
J11082570, J11084952, and J11112249. 

The qualitative spectral comparison in Fig. [5] shows that dust in the disks of cool stars 



-12- 



is more processed than it is in the ISM (see also lApai et al.ll2005l ). In fact, these disks 
have a broad, flat-topped 10 /im emission band, indicative of grains larger than the 0.1 jum 
grains that dominate the ISM. In addition, most spectra show peaks at 9.4 and 11.3 urn 



from crystalline silicates similar to those identified in spectra of comets (e.g., 



Wooden et al. 



20071 ). 



Van Boekel et al. (2003) first identified a correlation between the shape and the 
strength of the 10 /ttm silicate emission features from disks around intermediate-mass stars 
(Herbig Ae/Be stars). This correlation was interpreted in terms of dust evolution and 
is often used to characterize the amount of "primitive" versus "pr ocessed" dust in disk 
atmospheres. We have re-analyzed the Herbig Ae/Be sample from Ivan Boekel et i 
an d expanded the anal ysis to disks around young sun-like stars (from the sample studied 



20051 ) 



m 



Pascucci et al. 



20081 ) and cool stars (from this work). In Fig. 0, we show the flux ratio 



of the continuum-subtracted spectra at 11.3 and 9.8 /ttm versus the peak-over-continuum 
flux density within the 10 /im emission feature. Features dominated by "primitive" 
dust (sub-micron amorphous grains like those in the ISM) have low 11.3/9.8 ratios and 
high peak-over-continuum fluxes. In contrast, features dominated by "processed" dust 
(micron-sized grains and crystals) are characterized by high 11.3/9.8 ratio s and low 



p eak-o ver-cont inuum fluxes. This figure confirms the trend we reported in 



Apai et al 



( 120051 ): disks of cool stars have silicate emission features dominated by "processed" dust, 



the extent of pr ocessing is typically 



stars. Recently, 



l igher than in disks around sun-like and Herbig Ae/Be 



Kessler-Silacci et al. 



(120061 ) compared the strength of the silicate emission 
feature from disks around 13 Herbig Ae/Be stars to that from disks around 15 M stars 
and found a probability as large as ~15% that these two groups are drawn from the same 
parent population. We applied the Kolmogorov-Smirnov (K-S) test to our much larger 
sample of Herbig Ae/Be, sun-like and cool stars. We find that the feature strengths from 
the Herbig Ae/Be and sun-like star disks and from the sun-like star and cool star disks 



13 



are also consistent with being drawn from the same parent population (K-S probabilities 
of 6 and 7% respectively). On the other hand, the samples of cool star and Herbig Ae/Be 
disks statistically differ in the strength of their 10 /mi emission features (K-S probability of 
0.6%), with cool star disks having weaker features than Herbig Ae/Be disks. The shape of 
the 10 /xm emission feature is not statistically different among the three samples of disks 
(P(cool:sun-like)=7%, P(cool:Herbig)=16%, P(sun-like:Herbig)=6%). Note that the mean 
value of the peak-over-continuum becomes larger going from the disks of cool stars, through 
the sun-like stars, to the Herbig Ae/Be stars (see Fig. [6]). In Sect. 14.11 we will discuss a 
scenario that can explain the observed trend between the strength of the silicate features 
and the stellar mass/luminosity. 



3.2. Low HCN/C 2 H 2 Flux Ratios in Cool Star Disks 

The most exciting result of our survey has been the discovery of the vibration-rotation 
bands of C 2 H 2 and HCN in at least 5 out of 17 cool stars with disks (Figs. CEJ E], [3], |4j. 
No organic molecules have been previously detected toward such very low-mass stars. 
Motivated by the high detection rate, we have searched for the presence of the C 2 H 2 and 
HCN bands in the sun-like star sample. We find that at least 15 out of 44 sun-like spectra 
exhibit emission from C 2 H 2 and/or HCN (Figs. El [91 [10]). The sample of sun-lik e stars 



consists of about an equal number of single and binary stars (IPascucci et al 



20081 ). The 



detection rate of C 2 H 2 and HCN among single stars is similar to that among binary stars. 
Therefore, in the following we will not distinguish between single and binary stars, but 
rather contrast the sun-like star sample to the cool star sample. 

First, we note that the overall detection rate of C 2 H 2 and HCN is very similar in the 
sun-like star and cool star samples, with ~30% of the disks showing emission bands in the 
Spitzer low-resolution spectra. We searched for any correlation between the presence of gas 
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lines and the dust disk properties, such as the strength of the 10 am emission feature and the 

fl 

infrared slope of the source spectral energy distribution.]. Since the 44 sun-like stars were 
selected to have silicate emission features from optically thin disk atmospheres, we excluded 
the 3 featureless cool stars from this comparison (J11082570-7716396, Jl 1084952-7638443, 
Jill 12249-7745427). Fig. [TT] shows the histograms of the peak-over-continuum ratios of 
the 10/iin silicate emission features for the detections and for the non-detections of the 
C2H2 and HCN rovibrational bands. Detections are only present in disks with weak silicate 
emission features (peak-over-continuum ratios < 2). In addition, detections are confined 
to those disks with evidence for flatter geometries (Fig. [T2l) . In summary, the disks with 
gas detections show clear signs of grain growth and dust settling both among the sun-like 
and the cool star samples. This correlation may be an observational bias. Disks with larger 
grains and flatter disk geometries have we aker mid-infrared con tinuum levels making it 



easier to detect gas emission lines (see, e.g. 



Pascuccietal 



20071 ). 



When inspecting the detection rates of C2H2 and HCN separately, we find a striking 
difference between the samples of sun-like and cool stars. C2H2 is detected in at least 
5 out of 14 cool star spectra with silicate emission features, but at most in 7 out of 44 
sun-like stars (only 4 are considered 'firm' detections). In contrast, only the single cool star 
spectrum of Cha Ha2 presents a possible HCN emission band, but at least 13 sun-like stars 
have HCN detected (see Table |3]). In addition, the band strengths of C2H2 and HCN are 
different in the two samples. Emission from HCN is brighter than that from C2H2 in most 
disks around sun-like stars (a notable exception is CY Tau). In contrast, all disks around 
cool stars have brighter C2H2 emission compared to HCN. This trend is well illustrated in 
Fig. [HI where we plot the median (at each wavelength) of the continuum-subtracted and 



1 for this w e used the IRS fluxes at 7 and 13 yum, and the MIPS 24 yum data from 



Pascucci et al. 



(j2008|) and 



Luhman et al. 



(12008 
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normalized spectra for the sun-like and the cool stars with C2H2 and/or HCN emission band 
detections. The flux ratio of HCN over C2H2 is ~2.9 for the median sun-like star spectrum, 
while it is only 0.26 for the median cool star spectrum. The only two sources where both 
C2H2 and HCN are firmly detected are the sun-like stars DF Tau and CoKu Tau 3 with 
HCN over C2H2 flux ratios of ~1.3 and 0.7, respectively (see Fig. [14]) . These two objects 
alone point to a diversity in the line flux ratios of HCN and C2H2 in sun-like stars. This 
diversity should be further investigated with sensitive high-resolution Spitzer spectra. We 
stress, that the differences in detection frequency and band strength between the sun-like 
and cool stars are not explained by possible differences in the infrared continuum, its shape, 
and signal-to-noise ratio among the two samples because the C2H2 and HCN bands are at 
adjacent wavelengths. Sect. 14.21 suggests a possible scenario to explain these findings. 



4. Discussion 

4.1. Dust Processing, Stellar Luminosity and Disk Turbulence 

The main result of Sect. 13.11 is the trend between the stellar mass/luminosity and the 
strength of the 10 /mi silicate emission feature with disks around Herbig Ae/Be stars having 
statistically stronger features than disks around cooler/lower-mass stars. In other words, 
the 10 /an emission features from disks around low-mass/luminosity stars are dominated by 
more "processed" dust than the same features from disks around higher mass/luminosity 
stars. How can we explain this observational trend? 

The 10 jum silicate e mission feature arises from different regions in disks around stars 



of different luminosities. iKessler-Silacci et al. 



(120071 ) calculated that it probes radii > 



0.5-50 AU in disks around Herbig Ae/Be stars, but only radii <0. 001-0. 1AU in brown 
dwarf disks. This fact may provide an explanation to our observations. The closer the 
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grains are to the central star the faster they are depleted due to the higher density and 
rotational frequency which enhance the collision rate of dust grains and the velocity at 
which particles settle to the disk midplane. In this picture, the grain sizes inferred from 
the 10 /im silicate features are not representative for the bulk silicate dust in the disk but 
rather reflect the diff erent location of the 10 um emission zone a round stars of different 



luminosities (see also lApai et al. 



2005 



Kessler-Silacci et al. 



20071 ). However, grain growth 



and dust settling are not the only mechanisms determining the grain size distribution in the 
disk atmosphere. 

Turbulence may also affect the strength of the 10 /im emission feature, for instance 
by stirring up lar ge grains from the disk m idplane. In the frame of the magnetorotational 



instability model 



disk turbulence (IGammie 



Balbus fc Hawlev 



1996^ 1. 



199ll ). stellar accretion rates are a good proxy for the 



In Fig. [12] we plot the strength of the 10 /im emission 



features versus the stellar accretion rates for the sample of sun-like stars. The upper 
envelope of this plot suggests a trend with high accretion rate systems having large 
grains and moderate/low- accretion rate systems having small grains. However, the large 
uncertainty in the stellar accretion rates precludes any firm conclusion. This, in addition 
to the fact that large grains are seen commonly at all accretion rates, suggests that 
accretion/turbulence can only play a minor role in determining the strength of the 10 /im 
silicate emission feature. 

Finally, it should be also noted that disk models including only dust coagulation and 
settling/vertical mixing produce a t oo rapid (only 10 4 years) deple tion of grains up to 100 fim 



in size essentially at all disk radii (jDullemond fc Dominik 



20051 ). which is inconsistent 



with the frequent detection of 10 fim emission features in Myr-old disks. Fragmentation 
of aggr egate grains could b e the major mechanism replenishing the population of small 



grains. 



Brauer et al. 



(120081 ) implemented aggregate fragmentation in their disk model and 
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show that grain growth and fragmentation can reach an equilibrium for a period of ~ 10 6 
which results in a quasi- stationary grain size distribution. If the average grain size of 
quasi- stationary grain size distributions decreases with radial distance from the central star 
then such disk models could explain our observations. 



4.2. HCN under-abundance in Disk Atmospheres of Cool Stars 

We have shown in Sect. 13.21 that there is a striking difference in the detection rates 
of HCN and C2H2 between the sun-like and the cool star samples. HCN is often present 
in sun-like star spectra, but absent from the spectra of cool stars. Similarly, Fig. [141 
demonstrates that the HCN/C2H2 flux ratio is much higher toward sun-like stars than 
toward cool stars. As explained below, these differences arise from different abundances in 
the two samples. 

To model the emission from the v<s band of CoH? and the u 2 band of HCN we followed 



the procedure described in lLahuis fc van Dishoeck! (120001 ). which assumes rotational and 



vibrational levels in local thermal equilibrium. Using the populations in each level, the 
optical depths are calculated assuming a Voigt profile function. The spectrum is also 
multiplied at each frequency by the black body emission of gas at the temperature T gas and 
reduced to the resolution of the IRS/SL module using a Gaussian profile. First, we explored 
whether our low-resolution spectra can provide any interesting constraints on T gas . To do 
that, we have simulated spectra with gas temperatures from 100 to 900 K and measured the 
shift in the band head peak emission. The peak emission of the v<i band of HCN shifts from 
14.02 /j,m for T gas =100K to 13.96 //m for T gas =900K. This shift is about half the spectral 
resolution of the IRS/SL 1 module at 14 fim and cannot be reliably detected in our spectra. 
A similar result is reached for the z/ 5 band of C2H2. We find that the wavelengths at which 
the emission of HCN and C2H2 peak in the observed spectra are in between the predicted 
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peaks for gas at 100 K and 900 K. In summary, the resolution of our spectra is not sufficient 
to measure the temperature of the emitting gas, but it is fully consistent with the emitting 
gas having temperature between 100 K and 900 K. 

More excitingly, the HCN/C2H2 flux ratio informs us on the relative column densities 
of the two molecules. We computed synthetic spectra with varying column density ratios 
of HCN over C2H2, calculated the line flux ratios and compared them to the observed 
line flux ratios in the sun-like and cool star samples (see Fig. [T5l) . Absorption spectra of 
HCN and C2H2 toward young embedded massive stars find highly correlated excitation 
temperatures for these molecules demonstra ting that they probe warm gas at similar 



temperatures (ILahuis fc van Dishoeck 



20001 ) . We consider her e three temperatures for 



Carr fc Najita 



(120081 ) find to fit best the 



the HCN and C2H2 emitting gas: 650 K, which 
HCN and C2H2 band profiles from AA Tau, 450 and 850 K as comparison temperatures. 
Other important model parameters are the gas column densities and Doppler broaden 
line widths. In our calculations we keep the C2H2 column density fixed to 10 16 cm -2 
while we vary the H CN column density b etween 0.5-11 16 cm -2 , thus covering the column 



densities derived by 



Carr fc Najita 



(120081 ) for AA Tau. We also assume a Doppler line 
broadening of 5km/s, greater than the thermal broadening for the considered temperature 
range. The strongest C2H2 and HCN emission lines start to become optically thick around 
10 16 cm~ 2 and 7xl0 16 cm" 2 , respectively, for gas temperatures >400K (see Fig. 3 from 



Lahuis &: van Dishoeck 



2000j). Fig. [15] shows that, regardless of the gas temperature, 



HCN/C2H2 line flux ratios as high as 3 can be obtained only for large HCN/C2H2 column 
density ratios. The line ratio is almost independent from the gas temperature because both 
the C2H2 and the HCN lines start to become only marginally thick for the parameters 
we considered. If we assumed higher column densities for the two molecules or smaller 
Doppler line broadening, the HCN over C2H2 line flux ratio would start to saturate faster 
and even higher HCN/C2H2 column density ratios would be required to obtain the typical 
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HCN/C2H2 line ratio of 3 of sun-like stars. In contrast, the low HCN/C2H2 flux ratio 



measured in cool stars must originate from a 



C2H2. For comparison, 



Lahuis fc van Dishoeckl (120001 ) and 



o wer c olum n density of HCN in re spect to 



Boonman et al 



(120031 ) find that 



HCN is a factor of a few more abundant than C2H2 in hot molecular cores. How can we 
explain the HCN under-abundance in the disk atmospheres of cool stars? 

Here we propose a possible explanation. If HCN is formed in the warm disk 
atmo sphere then its abun dance may be limited by the availability of gas-phase atomic N 



e.g. 



Agundez et al. 



20081 ) . Atomic N is primarily produced via the photo-dissociation 



of N 2 , th ought to be 



; he m ost abundant nitrogen-bearing molecule in the interstellar 



medium (IBergin et al. 



19951 ). Thus, the HCN abundance may be limited by $ uv , the 



flux of UV-photons with A < 110 nm, required to dissociate N2LJ. For y oung stars $n V is 



dominated by emission o 



2003 



Herczeg et al. 



hot (~10,000K) gas accreting onto the stars (IMatsuyama et al. 
20041 ). A simple estimate for the accretion luminosities of the two 
samples suggests that their $ uv differ by more than 3 orders of magnitude, although the 
stellar temperatures differ by only 30% (Table [3]). If the HCN production from N is indeed 
driven by the UV-photochemistry, the large difference in the UV flux of sun-like and cool 
stars may naturally explain our observations. An important consequence of this scenario is 
the general under-abundance of molecular species produced through UV chemistry and an 
over-abundance of molecules wi th large bin ding energies in cool star disks. With efficient 



vertical and radial mixing (e.g., 



Ciesla 



20071 ). the composition of the disk atmosphere may 



be representative of the overall disk composition. Indeed, turbulent mixing seems to be 
necessary to explain the high abundances of warm HCN and C2H2 in the only single disk 



2 While the dissociation of NH 3 could also provide free N, Ivan Dishoeck et al.1 (120081 ) 
show that photodissociation rates for this molecule are insensitive to the radiation field for 
temperatures between 10,000 and 4,000 K 
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atmosphere where they were detected previously (ICarr fe Najitall2008l ). Detailed dust and 
gas disk modeling is required to understand whether the UV-photochemistry in the disk 
atmosphere can radically affect the organics in the disk midplane. 

An alternative, though perhaps less likely explanation is that our observations are 
tracing cometary or interstellar ice evaporation, in which case the organics in the disk 



midplane may not be different. 



DiSanti Mummal (120081 ) have shown that Oort cloud 



comets have an HCN/C2H2 ratio near unity, but Jupiter family comets have less C2H2. 
Because the HCN abundance relative to water is high by interstellar standards in Oort 
cloud comets, HCN seem to f orm efficiently in the outer disk, perhaps via surface chemical 



reactions on icy dust grains (jCharnley et al. 



1992l ) . Nonetheless, the evaporation of 



cometary ices alone cannot reproduce our results, unless C2H2 is preferentially destroyed 
in the disk atmospheres of sun-like stars. C2H2 has a dissociation cross-section that is at 
l east an o rder of magnitude higher than HCN for photons with A between 1300-2000 A 



flLee 



19841 ). Because sun-like stars have more flux than cool stars at these wavelengths, 
C2H2 could be more easily destroyed in the disk atmosphere of sun-like stars. However, it 
is interesting to consider how phqtodiss ociation rates depend on the stellar radiation field 



(Table 2 from 



van Dishoeck et al. 



20081 ). Going from a stellar radiation field produced by 
a star of 10,000 K to that of a 4,000 K star the C2H2 photodissociation rate decreases only 
by a factor of 9, that of HCN by a factor of 68, while that of N 2 by more than 4 orders of 
magnitudes. Based on this dependence on the temperature of the radiation field, it seems 
more likely that the large difference in HCN/C2H2 abundance is caused by more efficient 
dissociation of N2 in disk atmospheres of accreting sun-like stars. 
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4.3. On the Organics Available During Planet Formation 



HCN is the key building block of many complex organic compoun ds. Under reducing; 
conditions the polimerization of five HCN molecules can form adenine (IOrolll960l ). one of 
the two purine nucleobases at the core of the replication system of all terrestrial organisms. 
But even if the early Earth atmosphere was reducing, the synthesis of adenine thr ough HCN 



1993 ) and its 



polymerization under terrestrial conditions remains controversial (jShapirol 
delivery through come tary or meteoritic impacts has been proposed as a likely alternative 



( Chvba fe Sagan 



( IStoks fe Schwartz 



1992). Adenine synthesis 



1981 



Kissel fe Krueger 



abun dant HCN in the interstellar space 
2002) and recently in disks (this work 



likely occurs in extraterrestrial environments 



1987), consistent with the observations of 



Hirota et al. 



Lahuis et al.l 



2006 ; 



1998 ), in comets 



Gibb et al. 



(Magcc-Sau er et al. 



2007 



Najita et al 



20071 ). In addition, the availability of simple organic molecules such as HCN has a 



20041 ) 



fundamental impact on the chemical networks in protoplanetry disks (ISemenov et al. 
and thus on the compounds available during planet formation. 

Our observations show that HCN is under-abundant in the disk atmosphere of cool 
stars. With efficient mixing in the disk, this finding may reflect an overall deficiency of 
HCN in disks around cool stars. In this case, it is possible that the organics in the planet 
forming regions of cool stars, as well as the organics possibly delivered post-formation, will 
differ from those around sun-like stars. If exogenous HCN has played a key role in the 
synthesis of prebiotic molecules on Earth as proposed, then prebiotic chemistry may unfold 
differently on planets around cool stars. 
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5. Summary 

Using Spitzer/IRS we compared the gas and dust properties of 17 disks around cool 
stars/brown dwarfs to 58 disks around sun-like and intermediate-mass stars. The cool 
stars have spectral types later than M5, the sun-like stars between Kl and M5, and the 
intermediate-mass stars have spectral types earlier than AO. As the first large comparison, 
we found important and surprising differences between these samples. Our results can be 
summarized as follows: 

1. We identify 14 cool star disks with 10 fj,m silicate emission features. Features are 
broad and weak suggesting that grains at least a few-microns in size are present in 
the disk atmospheres. Most of the cool star spectra also show prominent crystalline 
silicate features. 

2. The sample of disks around cool and intermediate-mass stars statistically differ in 
the strength of their 10 ^m emission features, with cool stars having much weaker 
continuum-subtracted features than intermediate-mass stars. The observed difference 
demonstrates that the dust populations at the observed disk radii are different in the 
two samples. 

3. We report the first detections of organic molecules from disks around brown dwarfs. 

4. We find a striking difference in the detection rates of HCN and C 2 H 2 between the 
sun-like and the cool star samples, with HCN often present in sun-like stars, but 
absent from the cool star spectra. Similarly, the HCN/C2H2 flux ratio of the median 
sun-like star spectrum is an order of magnitude higher than that of the median cool 
star spectrum. We interpret these as evidence for different abundances of the two 
molecules in the disk atmosphere of sun-like and cool stars. Efficient photodissociation 
of N 2 in the disk atmosphere of sun-like stars may drive the higher production of 
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HCN. 
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Fig. 1. — Left Panels: IRS low- resolution sp ectra for the cool star sample. The IRAC 5.8 



and 8 /xm fluxes from 



Luhman et al 



(j2008al ) are over-plotted as filled black circles (when 
multiple observations are available we plot the mean flux values). Right Panels: Expanded 
view of the brown dwarf spectra around the C2H2 and HCN emission bands. A star indicates 



a "firm" detection while a question mark indicates a "possible" detection. 
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Fig. 4. — Same as Fig. [IJ 
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Fig. 5. — Qualitative comparison of the 10 /iin silicate emission features from cool star 
disks. Spectra have been continuum-subtracted and normalized t o the peak emission . For 



comparison the spectra of the amorph ous silicate-dominate d ISM (IKemper et al. 



the crystalline-rich comet Hale-Bopp ( ICrovisier et al. 



20051) and 



19971 ) are also shown. 
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Fig. 6. — Shapes and strengths of the 10 /im emission features for the samples of Herbig 
Ae/Be stars (green stars), sun-like/T Tauri stars (blue squares), and cool stars/brown dwarfs 
(red circles). On the y-axis we plot the ratios of the normalized flux at 11.3 /im over 9.8 fim, 
which are proxies for the degree of crystallinity. The x-axis gives the peak over continuum 
in the 10 /im region of normalized spectra, which indicates the amount of grain growth. The 
filled symbols with errorbars show the mean values and the la standard deviations of the 
peak-over-continuum and of the F11.3/F9.8 for the three groups of disks. Cool star disks 
have statistically weaker 10 /im features than Herbig Ae/Be disks. 
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The full spectra are presented in 



Pascucci et al. 



(120081 ) 
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Fig. 8. — Same as Fig. [7J 
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Fig. 10. — Same as Fig. [3 
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Fig. 11. — Upper Panel: Histogram of the 10/zm silicate feature strengths for the detections 
(solid line) and non-detections (dashed line) of the C 2 H 2 emission band in the combined 
sample of disks around sun-like and cool stars. Lower Panel: Histogram of the lOyum feature 
strengths for the detections (solid line) and non-detections (dashed line) of the HCN emission 
band in the sample of disks around sun-like stars (no cool star spectrum presents a firm 
HCN detection). Both panels show that detections are confined to weak features (peak- 
over-continuum ratios < 2). The K-S probabilities that C 2 H 2 and HCN detections and non- 
detections are drawn from the same parent population are only 0.3 and 2.8% respectively. 
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Fig. 12. — Upper Panel : Flux ratios at 24 over 7/im versus the signal-to-noise in the 
conti nuum close to the rov ibrat ional emission b ands o f C2H2 and HCN. Fluxes at 24 //m are 



from 



Pascucci et al. 



(120081 ) and 



Luhman et al. 



(j2008al ). Flatter disk structures are indicated 



by lower flux ratios in the figure. Detections of the C2H2 band are filled symbols, non- 
detections are empty symbols ("possible" detections are omitted). Squares are for sun-like 
stars, circles for cool stars. Lower Panel : Same as upper panel but for the HCN emission 
band. No cool star spectrum presents a "firm" HCN detection. C2H2 as well as HCN 
detections are confined to low flux ratios, indicative of flatter disk structures. The figure 
demonstrates that the pattern of the gas line detections are not primarily set by the signal- 
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Fig. 13. — Stellar accretion rate (M accr ) versus 10 feature strength 



sun- 



1998 



te stars. Accretion rates are from the literature (Harti gan et al. 



Gullbring et al. 



1998 



White k Basri 



2003 



a consistent scale using the factors calculated in 



Naiita et al. 


2007) 


Najita et al. 


(2007) 



1995 



or the sample of 



Hartmann et al. 



and have been placed on 



(120071 1 FM Tau, CoKu Tau 3, 



IT Tau, and V710 Tau do not have m e asure d accretion rates. The errorbar on the peak- 



over-continu um is from 
rate is from iNajita et al. 



Pascucci et al. 



(120081 ) while the uncertainty in the mass accretion 



(120071 ). There is a 4% probability that the data are randomly 
distributed, the negative Kendall rank correlation coefficient (-0.22) indicates a moderate 
anti-correlation. Yet, given the large uncertainty in the mass accretion rate we are cautious 
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Fig. 14. — Median of continuum-subtracted and normalized spectra for the sun-like star 
(blue) and the cool star (black) samples presenting C2H2 and/or HCN emission bands (5 
cool stars and 15 sun- like stars). The spectra are normalized to the peak of emission and 
scaled to match the C2H2 emission in the two samples. The errorbars are the standard 
deviations of the normalized spectra. If cool stars had the same flux ratio of HCN versus 
C 2 H 2 as the sun-like stars do, HCN emission would have been easily detected toward them. 
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Fig. 15. — Line flux ratio of HCN over C2H2 versus column density ra tio of the two molecules 
Circle s are values predicted from sythetic spectra calculated following lLahuis fc van Dishoeck 



( 120001 ) . The C2H2 column density is kept fixed at 10 16 cm 2 while the HCN column density 



va ries from 0.5-1 
by lCarr fc Najita 



xlO cm , thus covering the C2H2 and HCN column densities reported 



(120081 ) for AA Tau. Filled circles are for a gas temperature of 650 K while 
empty circles are for gas at 450 and 850 K. Arrows indicate the line flux ratios measured from 
the median cool star and sun-like star spectra in Fig. 14. The dashed lines report the HCN 
over C2H2 flux ratios for the only two sources in our sample where both molecular lines are 
firmly detected. The figure illustrates that the large difference in the HCN/C2H2 flux ratio 
between the cool star and the sun-like star samples translates into very di fferent column 



Carr &: Najita 



densit y ratios of the two molecules. We note that in the case of AA Tau 

j — L 

(120081 ) find an HCN/C2H2 column density ratio of 8, similar to what we require for the 
median spectrum of the sun-like stars. 
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Table 1. Sample of cool stars in the Cha I star- forming region identified via infrared 

photometry. 



Source 


2MASS J a 


Sp. Type 


Fe.7fim 






Other name 








(mjy) 


(mjy) 


(mjy) 




Hn 2* 


11034764-7719563 


M5 


11.4±0.8 


3.7±1.0 


9.1 


Chal 425.ISO 32 


CHXR 15* 


11054300-7726517 


M5.25 


13.2±1.1 


- 


8.5 


Chal 449.ISO 65 


ESO-Hq559 


11062554-7633418 


M5.25 


- 


- 


11.4 




ISO 79 


11063945-7736052 


M5.25 


8.0±2.0 


9.9±1.5 


8.2 




CHSM 9484 


11071181-7625501 


M5.25 


- 


- 


2.4 


Chal 711 


Cha Ha 1 


11071668-7735532 


M7.75 


6.7±1.8 


11.8±1.5 


8.1 


ISO 95 


Cha Ha 9 


11071860-7732516 


M5.5 


11.4±0.9 


9.2±1.3 


9.4 


ISO 98 


Cha Ha 2 


11074245-7733593 


M5.25 


24.5±1.4 


23.4±1.2 


23.3 


ISO 111 


ISO 138 


11081850-7730408 


M6.5 


15.3±0.8 


17.6±1.3 


1.8 


Chal 410 


J11082570- 7716396 


11082570-7716396 


M8 






1.0 




ISO 147 


11082650-7715550 


M5.75 


4.3±1.1 




4.7 




Cha Ha 6 


11083952-7734166 


M5.75 


12.7±1.5 




11.8 


ISO 152 


J11084952-7638443 


11084952-7638443 


M8.75 






0.8 




T37 


11085090-7625135 


M5.25 


7.5±1.5 




11.1 


ISO 157 


ISO 165 


11085497-7632410 


M5.5 


12±1.8 




10.1 




ISO 217 


11095215-7639128 


M6.25 


10.3±2.5 




20.7 


Chal 726 


ISO 252 


11104141-7720480 


M6 


5.5±0.8 


6.8±1.2 


8.1 




Hn 13 


11105597-7645325 


M5.75 


31.4±1.8 


27.5±1.9 


44.4 


ISO 259 


J11112249-7745427 i ' 


11112249-7745427 


M8.25 











Note. — All spectral types are from 



Luhman 



(2004); 



Luhman et al 



I 2008al l and 



Luhman et al 



1 20051) . They 



are determ i ned fr om optical spectra and have typical precisions of 0. 25. Fluxes at 6.7 and 14.3 (im are from 



Pcrsi et al 



I 2000), while fluxes at 8 ^m are from 



Luhman et al 



(2008a 



) . The 8 fim IR AC fluxes are mean values 



of fluxes from multiple campaigns. 

a Thc 2MASS source name includes the J2000 sexagesimal, equatorial position in the form: hhmmssss+ddmmsss 



Cutri et al 



2003) 



*The IRS spectra are consistent with photospheric emission, see text. 
b This source has no 8 fim IRAC data but has excess emission at 5.8 fira 
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Table 2. Log of the observations. 



Source 


Date 


Peak-up 




Exposure Times 






mm/dd/year 


Mode 


SL1 


SL2 


LL1 


Hn 2 


03/12/2005 


PCRS 


60s x 20 


- 


30s x 36 


CHXR 15 


03/12/2005 


PCRS 


60s x 20 


- 


30s x 36 


ESO-Ha559 


08/02/2006 


PCRS 


60s x 40 


60s x 12 


- 


ISO 79 


03/13/2007 


PCRS 


60s x 40 


60s x 10 




CHSM 9484 


08/02/2006 


PCRS 


60s x 40 


60sxl2 




Cha Ha 1 


03/12/2005 


PCRS 


60s x 20 






Cha Ha 9 


03/12/2005 


PCRS 


60s x 20 




30s x 36 


Cha Ho 2 


03/12/2005 


PCRS 


60s x 14 






ISO 138 


03/12/2005 


PCRS 


60s x 14 






Jl 1082570-7716396 


03/11/2007 


PCRS 


60s x 40 


60sxl2 




ISO 147 


03/11/2007 


PCRS 


60s x 40 


60sxl2 




Cha Ho 6 


03/12/2005 


PCRS 


60s x 20 




30s x 36 


J11084952-7638443 


03/11/2007 


IRS 


60s x 40 


60sxl2 




T37 


07/25/2006 


PCRS 


60s x 40 


60s x 10 




ISO 165 


08/02/2006 


PCRS 


60s x 40 


60s x 12 




ISO 217 


03/12/2005 


IRS 


60s x 20 




30s x 36 


ISO 252 


03/11/2007 


PCRS 


60s x 40 


60s x 12 




Hn 13 


03/11/2007 


IRS 


60s x 40 


60s x 12 




J11112249-7745427 


08/02/2006 


PCRS 


60s x 40 


60s x 10 
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Table 3. Statistics on the firm detections of HCN and C2H2 from the sun-like and cool 
stars with disks. Because the sun-like stars were selected to have silicate emission features, 
we have excluded from the cool star sample the 3 featureless cool stars (J11082570, 
J11084952, and J11112249). Also shown in the Table are the mean stellar temperatures, 
masses, radii, and accretion rates, used to co m pute the accreti o n luminosities. Values for 



individual stars are 



Kenyon fc Hartmann 



rom 



(199 



2005) 



Naiita et al. 



feQOJ); 



Luhman 



(2007); 



Muzerolle et al 



(2007); 



Pascucci et al 



White k Basri 



(120051 ) . 



(120031 ) 



(2008); 



Natta et al. 



Sample 


C2H2 


HCN 


< T* > 


< Mi, > 


< R+ > 


< M > 




detections 


detections 


[K] 


[M ] 


[Rq] 


M /yr 


Sun-like stars 


4/44 


13/44 


3,900 


0.8 


2.1 


1.3xl0~ 8 


Cool stars 


5/14 


0/14 


3,000 


0.1 


0.7 
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